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We propose a geometric growth model which interpo- 
lates between one-dimensional linear graphs and small- 
world networks. The model undergoes a transition 
from large to small worlds. We study the topological 
characteristics by both theoretical predictions and nu- 
merical simulations, which are in good accordance with 
each other. Our geometrically growing model is a com- 
plementarity for the static WS model. 
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1. Introduction 

Many real-life systems display both a high de- 
gree of local clustering and the small-world ef- 
fect [1121314115] . Local clustering characterizes the 
tendency of groups of nodes to be all connected to 
each other, while the small- world effect describes 
the property that any two nodes in the system can 
be connected by relatively short paths. Networks 
with these two characteristics are called small- 
world networks. 

In the past few years, in order to describe real- 
life small-world networks, a number of models have 
been proposed. The first and the most widely- 
studied model is the simple and attractive small- 
world network model of Watts and Strogatz (WS 
model) [BJ, which triggered a sharp interest in 
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the studies of the different properties of small- 
world networks such as the WS model or its vari- 
ations. Many authors tried to find more rigor- 
ous analytical results on the properties of either 
the WS model or on its variants that were eas- 
ier to analyze or captured new aspects of small- 
worlds [71819110111112113114] . The WS model and 
its variants may provide valuable insight into some 
real-life networks showing how real- world systems 
are shaped. However, the small-world effect is 
much more general, it is still an active direction 
of research to explore other mechanisms producing 
small- world networks. 

In real systems, a series of microscopic events 
shape the network evolution, including addition or 
removal of a node and addition or removal of an 
edge [16117118119120] . In addition, the number of 
nodes in some real-life networks increases expo- 
nentially with time. The World Wide Web, for ex- 
ample, has been increasing in size exponentially 
from a few thousand nodes in the early 1990s to 
hundreds of millions today. To our knowledge, all 
previous models of small- world networks either are 
static (not growing) or have focused solely on ad- 
dition of nodes and edges one by one. Therefore, it 
is interesting to establish a exponentially growing 
network model and investigate the effect of more 
events, such as removal of edges, on the topological 
features. 

In this paper, we present a simple geometric 
growth network model controlled by a tunable pa- 
rameter q, where existing edges can be removed. 
As in the WS model, by tuning parameter q, the 
model undergoes a phase transition transition with 
increasing number of nodes from a "large-world" 
regime in which the average path length increases 
linearly with system size, to a "small-world" one 
in which it increases logarithmically. We study an- 
alytically and numerically the structural charac- 
teristics, all of which depend on the parameter q. 

Our model is probable not as much a model of 
a real-world system as an example of vast variety 
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Fig. 1. Scheme of the growing network for the case of 
q = 0.5, showing the first four steps of the iterative process. 
The dashed lines denote the removed edges. 



of structure within the class of networks defined 
by a degree distribution, but it is the first grow- 
ing model that exhibits the similar phenomena as 
the famous WS model. Thus, it may be helpful 
for establishing more realistic growing small-world 
models of real-life systems in future. 



Now we compute the number of nodes and edges 
in N t . We denote the number of newly added nodes 
and edges at step t by L v (t) and L e (t), respec- 
tively. Thus, initially (t — 0), we have L V (Q) = 2 
nodes and L e (0) = 1 edges in No. Let N c (t) de- 
note the total number of internode intervals along 
the semi- ring at step t, then N C (Q) = 1. By con- 
struction, we have L v (t) — N c (t — 1) for arbi- 
trary t > 1. Note that, when a new node is added 
to the network, an interval is destroyed and re- 
placed by two new intervals, hence we have the 
following relation: N c (t) = 2 N c (t — 1). Consider- 
ing the initial value N c (0) = 1, we can easily get 
N c (t) = L v (t+l) = 2*. On the other hand, the ad- 
dition of each new node leads to two new edges and 
one old edge removed with probability q, which 
follows that L e (t) = (2 - q) L v (t) = (2 - q) 2 t_1 . 
Therefore, the number of nodes Vt and the total 
of edges E t in N t is 



V t = Y,L v (t j ) = 2 t + l 



(1) 



and 



2. The network model 



Et= J2L e (t J ) = (2-q)2 t + q-l 1 (2) 



The network model is generated by decimating 
semi-ring, which is constructed in an iterative way 
as shown in Fig. [TJ We denote our network after 
t (t > 0) time steps by Nt. The network starts 
from an initial state (t = 0) of two nodes, which 
are distributed on both ends of a semi-ring and 
connected by one edge. For t > 1, JV* is obtained 
from N t -i as follows: for each existing internode 
interval along the semi-ring of N t -i, a new node 
is added and connected to its both end nodes; at 
the same time, for each existing internode interval, 
with probability q, we remove the edge linking its 
two end nodes, previously existing at step t — 1. 
The growth process is repeated until the network 
reaches the desired size. 

When q = 1, the network is reduced to the one- 
dimension linear chains. For q = 0, no edge are 
deleted, and the network grows deterministically, 
which allows one exactly calculate its topological 
properties. As we will show below that in the spe- 
cial of q = the network is small- world. Thus, 
varying q in the interval (0,1), we can study the 
crossover between the one-dimension regular linear 
graph and the small-world network. 



respectively. The average node degree is then 

2E t ti „ , 3-2q 

lk) t = — - = (4 - 2o) ; - 

W V t y ' 2* + 1 



(3) 



For large t and any q, it is small and approximately 
equal to A — 2q. Notice that many real- life networks 
are sparse in the sense that the number of edges 
in the network is much less than V t {V t — l)/2, the 
number of all possible edges [21314] . 



3. Structural properties 

Below we will find that the tunable parameter q 
controls all the relevant properties of the model, in- 
cluding degree distribution, clustering coefficient, 
and average path length. 

3.1. Degree distribution 

The degree distribution P(k) is defined as the 
probability that a randomly selected node has ex- 
actly k edges. For q = 1, all nodes, except the ini- 
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tial two nodes created at step 0, have the same 
number of connections 2, the network exhibits a 
completely homogeneous degree distribution. 

Next we focus the case < q < 1. Let ki(t) 
denote the degree of node i at step t. If node i 
is added to the network at step U (ti > 1), then, 
by construction, ki(ti) = 2. In each of the subse- 
quent time steps, there are two intervals with one 
at either side of i. Each of the two intervals will 
create a new node connected to i, and each edge 
connecting the end nodes of the intervals could be 
considered to be deleted with probability q. Then 
the degree fcj(i) of node i satisfies the relation 

h{t) = h{t - 1) + 2(1 - q). (4) 

It should be mentioned that Eq. (|3|) does not hold 
true for the two initial nodes created at step 0. 
But when the network becomes large, these few 
initial nodes have almost no effect on the struc- 
tural characteristics. Considering the initial condi- 
tion ki(ti) = 2, Eq. (UJ) becomes 

fci(t) = 2 + 2(1 -q)(t-ti). (5) 

The degree of each node can be obtained explicitly 
as in Eq. (J5]) , and we see that this degree increases 
at each iteration. So it is convenient to obtain the 
cumulative distribution [4] 

oo 

Pcum(k) = W ( 6 ) 

k'=k 

which is the probability that the degree is greater 
than or equal to k. For some networks whose de- 
gree distributions have exponential tails: P(k) ~ 
e _fc / K , cumulative distribution also has an expo- 
nential expression with the same exponent: 

oo oo 

This makes exponential distributions particularly 
easy to spot experimentally, by plotting the cor- 
responding cumulative distributions on scmiloga- 
rithmic scales. 

Using Eq. (©, we have P cum (k) = J2v=k p ( k ) = 
P(f <r = t-^ y ).Hence 




Fig. 2. Semilogarithmic graph of cumulative degree distri- 
bution of the network ATi7 for different q. All the networks 
have size 131073. The solid lines are the analytic calculation 
values given by Eq. 1(8) . 

Pcum{k) = y ] — 
t'=0 * 

~ 2* + 1 + ^ 2* + 1 

t'=\ 

= 2~ 5 n^rr (8) 

The cumulative distribution decays exponentially 
with k. Thus the resulting network is an exponen- 
tial network. Note that many small-world networks 
including the WS model belong to this class [3]. 

In Fig. [2 we report the simulation results of the 
degree distribution for several values of q. From 
Fig. [2l we can see that the degree distribution de- 
cays exponentially, in agreement with the analyt- 
ical results and supporting a relatively homoge- 
neous topology similar to most small-world net- 
works |9ll2ll3ll4j . 

3.2. Clustering coefficient 

By definition, clustering coefficient [B] Cj of a 
node i is the ratio of the total number of edges 
that actually exist between all ki its nearest neigh- 
bors and the number ki(ki — l)/2 of all possible 
edges between them, i.e. Ci — 2ei/[ki{ki — 1)]. The 
clustering coefficient (C) of the whole network is 
the average of all individual C[s. 

For the case of q = 1, the network is a one- 
dimensional chain, the clustering coefficient of an 
arbitrary node and their average value are both 
zero. 
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For the case of q = 0, using the connection rules, 
it is straightforward to calculate exactly the clus- 
tering coefficient of an arbitrary node and the av- 
erage value for the network. When a new node i 
joins the network, its degree fcj and e; is 2 and 1, 
respectively. Each subsequent addition of a link to 
that node increases both k t and by one. Thus, 
ej equals to k t — 1 for all nodes at all steps. So one 
can see that there is a one-to-one correspondence 
between the degree of a node and its clustering. 
For a node v of degree k, the exact expression for 
its clustering coefficient is %, which has been also 
been obtained in other models [12113114121] , This 
expression indicates that the local clustering scales 
as C(k) ~ k~ l . It is interesting to notice that a 
similar scaling has been observed empirically in 
several real-life networks [22] . 

Clearly, the number n(C, t) of nodes with clus- 
tering coefficient C = 1, |, |, • •• , j^j, \, , 
is equal to L v (t), L v (t — 1), L v (t — 2), • • • , L v (2), 
L v (l), L v (0), respectively. Therefore, the cluster- 
ing coefficient spectrum of nodes is discrete. Us- 
ing this discreteness, it is convenient to work with 
the cumulative distribution of clustering coeffi- 
cient [2H as 



W cum (C) = V n(C',t) = 2 

C'<C 



1 



c 



(9) 



where C and C are the points of the discrete spec- 
trum. The average clustering coefficient (C) t can 
be easily obtained for arbitrary t, 



(C)t 



1 

v t 
1 

1 • - 
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(10) 



For infinite t, (C) = - Ln(l - |) = In 2 , which 
approaches to a constant value 0.6931, and so the 
clustering is high. 

In the range < q < 1, it is difficult to derive an 
analytical expression for the clustering coefficient 
either for an arbitrary node or for their average. 
In order to obtain the result of the clustering coef- 




Fig. 3. The clustering 
a function of q. 



coefficient C of the whole network as 



ficient C of the whole network, we have performed 
extensive numerical simulations for the full range 
of q between and 1. Simulations were performed 
for network Nn with size 131073, averaging over 
20 network samples for each value of q. 

In Fig. [31 we plot the clustering coefficient C as a 
function of q. It is obvious that C decreases contin- 
uously with increasing q. As q increases from to 
1, C drops almost linearly from 0.6931 to 0. Note 
that although the clustering coefficient C changes 
linearly for all q, we will show below that in the 
large limit of q, the average path length changes 
exponentially as q . This is different from the phe- 
nomenon observed in the WS model where C re- 
mains practically unchanged in the process of the 
network transition to a small world. 

3.3. Average Path Length 

We represent all the shortest path lengths of N t 
as a matrix in which the entry dtj is the geodesic 
path from node i to node j, where geodesic path is 
one of the paths connecting two nodes with mini- 
mum length. A measure of the typical separation 
between two nodes in Nt is given by the average 
path length d tl also known as characteristic path 
length, defined as the mean of geodesic lengths 
over all couples of nodes: 



(11) 
(12) 



V t (V t - l)/2 ' 



where 
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denotes the sum of the chemical distances between 
two nodes over all pairs. For general q, there are 
some difficulties in obtaining a closed formula for 
dt- For q — 1 and q — 0, the networks are deter- 
ministic, which allows one to calculate dt analyti- 
cally. 

3.3.1. Case of q = 1 

For this particular case, the network is a linear 
chain (graph) which has two nodes with degree 1 
at both ends of the chain and Vt — 2 nodes with 
degree 2 in the middle. For convenience, we denote 
the total distances between all pair of nodes and 
average path length of a linear chain with n nodes 
as cr(n) and l(n), respectively. We label each node 
in the linear graph with size n from one end to the 
other as d = 1, 2, • • ■ , n — 1, n. Then we have the 
following equation: 

<j(n + 1) = a(n) + 9 n , (13) 



A 




Fig. 4. For q = 0, the network after t + 1 step iterations, 

N t +i, is composed of two copies of Nt denoted as and 

(2) 

iv t , which are connected to each other as above. 



where At is the sum over all shortest paths whose 
endpoints are not in the same Nt branch. The so- 
lution of Eq. {HI) is 



n(n + 1) 



where 9„ is defined as 

n n 

On = ^ ^i(n+l) = E 
i=l i=l 

Then the solution of Eq. {JT3J) is 

n— 1 n—1 

a(n) = a{2) + d rn = 1+13 



(14) 



m(m + 1) 



m=2 
'n-l 



l + \ E™ 2 + E 

\m=2 

n(n — l)(n + 1) 



m=2 

TO 
m=1 ) 



G 



Thus 



I(n) = 



a(n) 



n+l 



n(n-l)/2 3 ' 
which increases linearly with network size. 



(15) 
(16) 



3.3.2. Case of q = 

In the special case, the network has a self-similar 
structure that allows one to calculate d t analyti- 
cally, based on a similar method as that proposed 
in [52]. As shown in Fig. HI the network N t +i may 
be obtained by joining 2 copies of Nt , which are 
labeled as 7V 4 (q) , a = 1,2. Then we can write the 
sum St+i as 



S t+ i =2S t + A t , 



(17) 



St = 2*- 1 S'i + 2 t_m_1 A r 



(18) 



The paths that contribute to A t must all go 
through at least one of the 3 edge nodes (A, B, D) 
at which the two Nt branches are connected. Be- 
low we give the analytical expression for A t , called 
the crossing paths. 



We define 



d\ ot 



d „oar 



jmid 



E d%A > 

E diA , 



lEN. 



(i) 



diA<diD 

E * 

diA=diD 



y-near 



ymid 



jfar __ 



E diA > 



far 



E i. 

diA<d iD 

E l 
E 



d iA >d iD 



d iA >d iD 



so that c£ ot = d? ear + d™ id + c^ ar and V t = V f n 



rmid i T/far 



ymid + ytar^ By symmetry yneax = ^ ar< Xhus 
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A 



D 



Fig. 5. The first three steps of the network with q = 0, with 
the nodes labeled by ordered pairs denoting the shortest 
distance to the left and right end nodes. 



A, 



E 



^ [diA + d Aj 



E 



i,jj£A, d iA <diD 

{d iA + d A j) 



E 



(d lD + 1 + dsj) 



diA>dijji djA>djB 



[(v t -i)d lA + 4 ot ] 



ieJv t (1) , ijtA 

diA<d iD 



+ E [( v t acaz + v t mid - !)diA + 4 lcar + df d ] 



diA>dio 



]T [V t neaT (d iD + l) + <% eaT ] 



i£Nf>, i^A 
diA>diu 



dered pairs dnj, diA, for the first three iterative 
steps. We denote by a l m m , the number of nodes 
added at step t which have diA = m, diD = m! ■ 
Since A and D are connected, ml and m can dif- 
fer by at most 1. Thus for a given to there are 
three categories of nodes added at step t, respec- 
tively numbering a 



Having A t in terms of the quantities in Eq. (TT^) . 
the next step is to explicitly determine these quan- 
tities. We consider a node i € and the 
shortest-path distances to A and D, i.e., diA and 
dn). After the step when the node i was gen- 
erated, the values of diA and din do not change 
at subsequent steps, since the shortest paths are 
always along the edges added earliest. We see this 
in Fig. where the nodes are labeled by the or- 



m,m+b a m,m! an< ^ a m+l,m- By 



symmetry a* n+1 , m 



The to, to' values 



of nodes added at step t depend on the neighbor 
nodes, which were added at previous steps. For ex- 
ample, there is one node added at step t (t > 2) 
which is a nearest-neighbor of A, so this new node 
has to = 1, to' = 2, giving a\ 2 = 1. Nodes with 
to = 1, to' = 2 will in turn get neighbors with 
to = 2, to' = 3 in subsequent steps. The relation- 
ship between a\ 3 and a\ 2 for k < t is 



Similarly, 



t-2 

"2,3 = 

k=2 



t-2 



2 2 = 2(t - 3) . 



(21) 



t 3,4 



= E 2fl 2,3 = 2(<-4)(i-5). 



(22) 



Since nodes with distances to, to+ 1 do not appear 
before the step t = 2m, the sum over 3 starts at 
k = 4. Proceeding in this way, for general to > 1 
and i > 2to, 



= (F t - l)(d" car + d? id ) + (K" car + K mid - lK ot 
+ (F t ncar + V t mid - l)d l t ai + v t acai (d^ + df d ) 
+ ear (df ar + V t ncar ) + ^ ncar dJ iear . (20) 



l m,m+l ~ y ] 2a m _ l m 
fe=2(m-l) 
[>m-l 



(t — TO — 1)! 

(to - !)!(£ - 2m)! 



(23) 



The value of dp x is 1 for t = and for i > 0. 
Analogously, for general to > 2 and t > 2m — 1, 



L (t-TO-1)! 



(to-2)!(<-2to+1)! 



(24) 



The value of a\ 1 is 1 for t = 1 and for t > 1. 

As derived in the Appendix, we can obtain the 
quantities in Eq. [fTO]). 
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* Lfc/2J f lot, i Hsodd 

TAncar _ fc _ J 3 2 + 3' 1 1S 0Cta 

— °m,m+l — 1 1 

fc=l m=0 



6000 



3 2* + |, t is even 



v; mid = J2 E ak — J s 2 * + s> iisodd 

fc=l m=l 



t Lfe/2J 

EE 



ma 



m,m+l 



fe=l m— 



I ~ IT 2 * ~ ¥ + I* 2 *' * is odd 
- ^2* + it + \t2\ t is even 

t Lfe/2J 



^f d = E E 



k=l m=l 

__ im + i? 2 * + lt+h t2t > t is odd 

" \ii + ^2 i -ii+^2 t , t is even 

(25) 

where [^J denotes the largest integer < x and the 
different results for t odd and even are given con- 
secutively, and 



t Lfc/2J 

di r = E E ( w + iX^+i = d " car + ^ t ncar ■ 



fc=l m=0 



(26) 



Substituting the results of Eq. (J2SJ) into E q- CT - 
we can otain 

A t = — [-5 - 3(-l)* + 12 • 2* + 14 • 2 2t + 12*4*1 
18 

(27) 

Substituting Eq. J27D for A m into Eq. JTSJ), and 
using 5i = 3, we have 



f • 2* + i (4 - 5 • 2 2t - 10 • 2* + 6 1 2* 



-6t2 2 *),t is odd; 
+ 

+6t2 2t ),t is even. 



| ■ 2* + ^ (6 - 5 ■ 2 2i - 10 • 2* + 6 * 2* 



(28) 

Inserting Eq. (|28p into Eq. (jlip. we obtain the ex- 
act expressions for average path length dt which is 
of the form 

f 17+2 2 - t -5-2 t +6t(l+2 t ) / isnr ld- 

7 _ I 9(1+2^ . * 1S ° dd > f2ql 

' 1 17+3-2 1 -*-5-2 t +6t(l+2 t ) , . ^ J > 

1 -nr^r — i i, t is even. 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0. 

q 



9 1 



Fig. 6. Graph of the dependence of an average path length 
on the tunable parameter q. 



In the infinite network size limit (t — ► oo), 



-t 



\a.V t 



(30) 



9(1+2*) 



Thus, the average path length logarithmically 
grows with increasing size of the network. This log- 
arithmic scaling of d t with network size, together 
with the large clustering coefficient obtained in the 
preceding section, shows that in the case of q = 
the graph is a small- world network. 

3.3.3. Case of < q < 1 

For < q < 1, in order to obtain the varia- 
tion of the average path length with the parameter 
<?, we have performed extensive numerical simula- 
tions for different q between and 1. Simulations 
were performed for network N±$ with size 16385, 
averaging over 20 network samples for each value 
of q. In Fig. [51 we plot the average path length as 
a function of q. We observe that, when lessening q 
from 1 to 0, average path length drops drastically 
from a very high value a small one, which predicts 
that a phase transition from large-world to small- 
world occurs. This behavior is similar to that in 
the WS model. 

Why is the average path length low for small ql 
The explanation is as follows. The older nodes that 
had once been nearest neighbors along the semi- 
ring are pushed apart as new nodes are positioned 
in the interval between them. From Fig. [1] we can 
see that when new nodes enter into the networks, 
the original nodes are not near but, rather, have 
many newer nodes inserted between them. When 
q is small, the network growth creates enough 
"shortcuts" (i.e. long-range edges) attached to old 
nodes, which join remote nodes along the semi-ring 
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one another as in the WS model [6J. These short- 
cuts drastically reduces the average path length, 
leading to a small- world behavior. 



itive integer). For odd t, we can proceed in the 
same way. 

First, we compute y^ near . We rewrite m+1 and 



4. Conclusions 

In summary, we have presented and studied a 
one-parameter model of the time evolution of geo- 
metrically growing networks. In our model, in ad- 
dition to the widely considered case of addition 
of new nodes and new edges which connect new 
nodes and old ones, we also consider edge (linking 
old nodes) removal. The presented model interpo- 
lates between one-dimensional regular chains and 
small-world networks, which allow us to explore 
the crossover between the two limiting cases. We 
have obtained the analytical and numerical results 
for degree distribution and clustering coefficient, 
as well as the average path length, which are de- 
termined by the model parameter q. The observed 
topological behaviors are similar to those in the 
WS model. 

Our model may provide a useful tool to inves- 
tigate the influence of the clustering coefficient 
or average path length in different dynamics pro- 
cesses taking place on networks. In addition, us- 
ing the idea presented here, we can also establish 
a scale-free network model which undergoes two 
interesting phase transitions: from a large world 
to a small world, and from a fractal topology to a 
non-fractal structure. The details will be published 
elsewhere. 



Acknowledgements 

This research was supported by the National 
Natural Science Foundation of China under Grant 
Nos. 60373019, 60573183, and 90612007. 



„t _ nm-l/<ra-l n < 4-. 



_ 9 m-2r»m- 2 .! , , (31) 

a m,m ~~ z u * 



-t-m-U 2m<t + l 
For even t = 2T, we have 



= E E <m+i(m >1) + 1(A itself) 

fc=l 2m<k 

T 2T 

~ E E a m,™+l + 1 
m—l k—2ra 

T 2T 

= E E 2 " l ~ lc r;>-i + 1 

rri — 1 k—2m 
T 2T 

-E 2 ™" 1 E T-^-i + l- (32) 

m=l k=2m 

To find X)fe=2m ^feXi-l> we use tne approach 
based on generating functions [2"4"j . We define 

2T 

G 1 (x)= £ (l + x) k ~ m - \ (33) 

k—2m 

and denote by z the coefficient of power of £ m_1 
in G x {x). Then, we have z = Y%Li m T-m-r At 
the same time, we can sum the items in right hand 
of Eq. (|33|) and obtain G\ (x) as 



Gx(x) 



1/11 „,\2T-2m+l 

+ : (1 + x)- 1 

1 - (1 + x) 

(l + x) 2T - m - (1 + x) m - 1 



(34) 



Appendix 



From Eq. (|34|) . it is obvious that the coefficient 
of power of x m ~ x in G\(x) is C^_ m , thus z = 

^.CW! = C 2 ^_ m . Then Eq. © can be 
written as 



A. Derivations of V t neal and d% ear 

In this Appendix, we give the computation de- 
tails of y t noar and df car , respectively. Other quan- 
tities in Eq. (fl~9| can be derived analogously. Here 
we only address the case of even t — 2T (T is pos- 



v? 



^2™- 1 C 2 m T _ m + l. (35) 



All that is left to obtain V' t near is to evaluate the 



sum in Eq. ([35]), which is denoted as u. In order 
to find it, we define 
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G 2 (x) = J2 2 m " 1 (l + -' 



x .2T-m x T-m 



(36) 



Then u exactly equals the coefficient of the power 
of x T in G 2 {x). By summing over all m, we get 
G 2 (x) as 



G 2 (x) 



2 T (l + x) T -x T {l + x) 
(l-x)(x + 2) 

1/1 1 



2T 



3 \1 -x x + 2 / 
[2 T (l + x) T -x T (l + x) 2T ] 



And since 



1 - x 



= 1 + x + x z 



t=0 



and 



1 1 



a: + 2 2 ^ V 2^ 

t=o 



(37) 



Qt=Y, rn2 m - l C? T _ m 

m— 1 

= ^2— 1 (2T-m)CST 1 m _ 1 

m— 1 

= ^2— 1 (2T-l)C-" 1 m _ 1 

m— 1 
T 

-^2™- 1 (m-l)C 2 m T - 1 m _ 1 

= i(2T-l)(2 2T - 1 + l)-2Q T _ 1 . (40) 

Making using of Qi = 1, Eq. (|4"0")) can be solved 
inductively, 

Q T = (-2f- 1 Qi 

+ Y^^2) T - m - ll -{2m + 1) (2 2m+1 + 1) 
, 3 

m—l 

= — — 4 T + -T+-T4 T (411 

27 27 9 9 ' [ } 

as given in Eq. (f2"5|) . 
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t 
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T 

= £ m2 m - 1 C 2 m T _ m . (39) 
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